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Abstract 

For the system of Maxwell equations of electromagnetism in an Z-periodic composite medium of overall 
size L (0 < I < L < oo), in the low-frequency quasistatic approximation, we develop an electromagnetic 
version of strain-gradient theories, where the magnetic field is not a function of the magnetic induction 
alone but also of its spatial gradients, and the electric field depends not only on the displacement but 
also on displacement gradients. Following the work (Smyshlyaev, V.P., Cherednichenko, K.D., 2000. On 
rigorous derivation of strain gradient effects in the overall behaviour of periodic heterogeneous media, J. 
Mech. Phys. Solids 48, 1325-1357), we develop a combination of variational and asymptotic approaches 
to the multiscale analysis of the Maxwell system. We provide rigorous convergence estimates of higher 
order of smallness with respect to the inverse of the “scale separation parameter” L/l. Using a special 
“ensemble averaging” procedure for a family of periodic problems, we derive an infinite-order version of 
the standard homogenised operator of second order. 


1 Introduction 

Heterogeneous media with microstructure (“composites”), which in the simplest case involve two length- 
scales: “macroscopic” L and “microscopic” I, so that L ^ I, have been a focus of attention across engi¬ 
neering and physical sciences since the advent of quantum mechanics and the subsequent development of 
the ideas of “upscaling” and “course-graining” in the quantitative description of material properties. It is 
widely known that under the assumptions of periodicity (where both L, I are periods in the variation of 
material properties) and strict separation of scales, i.e. in the regime when Z/L —>■ 0 is a suitable approx¬ 
imation, the behaviour of composites is fully determined by the solution of the so-called “cell problems” 
on the period of the composite. 

It has also been noticed in the context of elastic solids (see e.g. m) that when the ratio l/L is “small 
but not too small”, additional terms need to be taken into account in order to obtain an accurate picture of 
the overall response of the composite. Mathematically speaking, the solution to the cell problem does not 
suffice when the length-scales involved are not well separated from each other, and a more sophisticated 
“higher-order” averaging framework has to be invoked, as was shown in m, 0. 

Keywords: Homogenisation, Maxwell equations, constitutive laws, asymptotic analysis. 
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In the present work, we develop the higher-order approach for the analysis of the classical system of 
Maxwell equations. Our motivation for doing so stems from the increasing interest to “metamaterials” 
(see e.g. [IS] for an overview) and from the indication (see m. 0) that the quantitative description of 
the effect of length-scale interactions on the overall behaviour of composites is key to developing new ways 
of manufacturing metamaterials. 

The centrepiece of the homogenisation theory for second-order elliptic equations (see e.g. m, m, 
[I]) is the derivation of the homogenised equation that captures the “effective” behaviour of the original 
problem by “averaging out” small-scale oscillations. By contrast, the main products of our approach are 
the “infinite-order homogenised equation”, see ()3.19p . (I4.21h . involving a sum of higher-order gradients 
of the classical homogenised solution, with coefficients given by increasing powers of the scale parameter 
£ := l/L, from which the classical homogenised equation is obtained as the truncation of lowest order, and 
the “homogenised equation of higher order” obtained as the Euler-Lagrange equation for the “higher-order 
homogenised variational problem” (I4.13p . 

The use of asymptotic expansions in understanding size effects in periodic media was proposed in |3|, 
m, where the solution to linearised elasticity equations in a periodic composite medium was sought in the 
form of a two-scale asymptotic series expansion whose terms depend on the macroscopic x and microscopic 
x/e variable. Using a special averaging procedure, a set of “higher-order stress-strain relations” for such 
media is derived in |3] and m, where the stress is a function of not only strain but also of gradients of 
strain. The advantage of the asymptotic approach in deriving such constitutive laws is the availability 
of rigorous error estimates, for small values of e, for the difference between the actual strain and stress 
and their higher-order versions. However, this approach also has two drawbacks from the perspective of 
numerical implementation the higher-order expansions are only expected to be accurate if e is sufficiently 
small and may lead non-elliptic differential equations, as shown in |1|. 

Following [12], who study the equations of elasticity, we develop a higher-order homogenisation frame¬ 
work for the system of Maxwell equations, via a combination of asymptotic and variational calculi. In the 
first part of our analysis. Section [3l we provide a higher-order extension of the classical two-scale asymptotic 
approach, including the higher-order version of the standard convergence estimates. In the second part 
(Section Hj), we use a version of “variational asymptotics” to derive higher-order homogenised equations 
that are elliptic by construction. This is achieved by using a set of trial fields suggested by the asymptotic 
analysis in a family of variational formulations associated with the problem for the original composite 
medium. Furthermore, we show that the corresponding higher-order variational solution is “close”, in a 
certain variational sense (Proposition 14. 3|) . to the solution of the original problem. The main idea, which 
we adopt from m, is to cancel the effect of the rapid oscillations in higher-order terms via the “ensemble 
averaging” of a family of problems obtained by shifting the fast variable of the original problem. 

By analogy with m, the asymptotic and variational approaches result in equivalent “infinite-order” 
homogenised equation (Sections 13.11 and 14.4p . however the proof of this fact in the Maxwell case requires a 
special tensor symmetrisation procedure, which we carry out in Section [5| 

Finally, in Section [6l we derive the infinite-order effective constitutive relations between the magnetic 
field and induction and between electric field and displacement. In analogy with the elasticity case, the 
magnetic field depends not only on the magnetic induction but also on its spatial derivatives, with coeffi¬ 
cients that become more significant for larger values of e. This leads to the generalisation of the notion of 
macroscopic magnetic permeability, in the form of an “effective permeability operator”, whose inverse is a 
differential operator of infinite order, see Section [6.31 

In the main body of our work (Sections [2H4]) we focus on the pair of equations governing the behaviour 
of the electric field when the electric permittivity is constant. Similar analysis, however, applies to the pair 
of equations for the magnetic field, where electric permittivity rapidly oscillates and magnetic permeability 
is constant, subject to additional considerations due to the presence of an oscillatory current density term, 
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as discussed in Section 16.41 By analogy with the case of the magnetic field and induction, we derive an 
infinite-order “effective permittivity operator”, which links the overall electric field and displacement. 

2 Formulation of the problem 

In what follows we study the vector equation 

curlj^^—^curl u'^(x)| = f(x), x€T:=[0,T]^, e,T>0, T/eGN, (2.1) 

where the coefficient matrix A is assumed to be measurable, Q-periodic, symmetric: 

Aij(y) = Aji(y) Vy€Q:=[0,l)^ i,j = 1,2,3, 
bounded and uniformly elliptic: 

3i/>0: < Aij{y)^i^j < V^GM^y€Q. 

Throughout the paper, we use the notation A^{-) := A(-/e) for matrix functions A with the above proper¬ 
ties. We take right-hand sides f that are infinitely smooth in M^, T-periodic, divergence-free and have zero 
mean. The equation ()2.1I) describes the behaviour of the electric component of an electromagnetic field in 
the quasistatic approximation (see Section [6^ . Here the quantities A, W, and f represent the inverse of 
the magnetic permeability /I, the electric field Ef and the current density —Jo, respectively, at each point 
X G M^(modT). In terms of the “size parameters” l,L mentioned in the previous section, we assume that 
e/T = l/L. 

Define the space to be the closure of the set [C'j^j.(T)]^ of infinitely differentiable T-periodic 

functions with respect to the norm 



|curlu||[i2(T)]3. 

(2.2) 

We shall study weak solutions W to (|2.1l) in the space 



T(T) := {u G [t 2(T)]^ divu = 0, 

(u)t = 0} n FeurlW, 

(2.3) 

so that the identity 

r r 



/ T'^curl W • curl (^ = / 1 ■ Lp 

Jf Jt 

V(^G [C-p-(T)]' 



holds. The divergence div in (12.31) is understood in the sense of distributions. Henceforth, we use the 
notation ( • )t, ( ' ) for averages (mean values) over T and Q, respectively: 


and the averages of vector quantities are taken component-wise. Note that T’(T) equipped with the norm 
()2.2I) is a Sobolev space. For all e > 0, the equation (|2.ip is well posed in X(T), as shown next. 

Theorem 2.1. For all f € ^ hi T’(T), there exists a unique solution W G T’(T) to 12. 1\} . 
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Proof. Define a bilinear form b on ^(T) by the formula 


b(u, (^)= / A*^curl u • curl u, 

■Jt 

A unique solution to the problem 

b(u^,(^)= [i-p 'iip ^ Xiff) 

Jt 

exists by the Lax-Milgram Lemma (see e.g. ID)- Indeed, coercivity of the form b follows from ellipticity 
of the matrix A : 


b(u, u)> f z^jcurlu |2 > iz. / 

|curlup + —z/ f |curlup > —u f |curlup + -z/C'i f |up > C 2 ||u| 
Jt 2 Jtt 2 Jt 2 Jt 2 Jt 


^curl(T) 


for some Ci,C 2 > 0, where we have used (16.251) . see Appendix A, combined with the fact that u is 
divergence-free. Continuity of the form b follows from the Cauchy-Schwarz inequality and the assnmption 
that A is bounded: 


b(u,(^) < ( / IA curl u 


1/2 


I curl (p\ 


9\V2 - , „ 

' < Cl' ^ U TTl 




for some C > 0. 


□ 


3 Asymptotic expansion of the solution to (l2.1l) 


We seek a solution of ()2.ip in the form of a two-scale power series: 


= ^e^Uj(x, ^), xGT. 
j=o ^ 

Substituting (13.ip into equation (j2.ip yields a more specific form for the coefficients Uj, namely: 


(3.1) 


W(x) = v(x,e)-F^e'^|Vy(^(^)(y)V4v(x,e))-h ^)(y)V4 ^v(x,e)) 


i=i 




-1 


curlxv(x, e)| 


y=x/£ 


(3.2) 


where we set = 0. Henceforth, we denote by Vx, Vy, j = 0,1,2,..., the tensors of derivatives of 

order j with respect to the variables x, y, often omitting the lower index when it is clear from the context. 
The coefficients and are tensors of order J -|- 1, J = 0,1,... , whose components belong to the 

spaces iJpgj. g(Q) and respectively, where iLpgj. Q((5), Z = 1,2, denotes JJ*((5)-closure of the set 

of elements of [^“^.(Q)]^ with zero mean. The smooth, T-periodic, divergence-free vector field v(-,e) is 
sought as a series in powers of s : 


v(x, ^)=Y 1 W’ X G T, 
k=0 


(3.3) 


where v^. A: = 0,1, 2,... , are e-independent. 
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Remark 1. The tensor products in (13.21) are evaluated as follows: 






^curlxV = ^S,„i.^,(curlxv)i.^,,i2...i., 


for ifc G {1, 2, 3}, k = 1,2,... , j + 1, where summation is carried out for repeated indices, and the comma 
denotes differentiation with respect to the indices following the comma. The curl and divergence of a tensor 
of order j + 1 are tensors of order j + 1 and j, respectively, given by 

(curloT'*^-^^). . . = id\Y . . . = ,■ o- 

We formally substitute (|3.2p into the equation (j2.ip . treating the “slow” variables {xi,X 2 ,xz) =: x and 
the “fast” variables ( 2 / 1 , 2 / 2 ) ys) =: y independently, so the “full” gradient and curl operators are evaluated 
according to the rules V = Vx + e“^Vy, and curl = curlx + curly. Making use of the identities 
curlxVx(') = curlyVy(-) = 0 and curlxVy(-) = —curlyVx(-): we formally write 

00 

curlx{^ curlxv} + e“^curly curlxv} + ^e^curlx{Acurlx(oT'^^^V4 ^curlxv)} 

3=^ 

00 

+ ^curlyj Acurlx(oT'^-^^V^“^curlxv) | + curlx{M curly (M^^-^^V^“^curlxv) 

00 

+ g'^~^curly |M curly (c/h^'^^Vj^~^curlxv) | = f. (3.4) 
i=i 

Equating the terms in p3.4p corresponding to individual powers of e, we obtain: 

curlyjAcurly (M^^^^curlxv) I = —curly {^4 cur lx v}, (3.5) 

curly {j 4 curly (M^*^^^VxCur lx v)} = f — curlx{^ curlxv} 

— curly {M curlx curlx v)} — curlx curly cur lx v)}. (3.6) 

Requiring that (j3.5p be satished for all admissible vector fields v implies 

curljAcurlM^^^^} = —curl A, (3.7) 

which is understood in the weak sense, i.e. 

((AcurlM^*^^) + A)curl(/)) = 0 V(/) G [C'“j.((3)] (3.8) 

The matrix is determined uniquely under the conditions that it is Q-periodic and has zero average. 

The condition of solvability of (13.6p . viewed as an equation for is the usual homogenised equation 

(c/. [I]) 

curlxjh^^^curlxv} = f, vGT’(T), := (y4(curl +/)). (3.9) 

Remark 2. In the case of a “laminate”, ^(y) = a{y 2 )I, the matrix is found to have the form 


0 0 -iV^ 

=(00 0 

WOO 


(3.10) 


where N = N{y 2 ) satisfies —{aN')' = a', subject to the conditions that it is (^-periodic and has zero mean. 
Direct calculation shows that 

r-l |^y2 

((a“^)“^Q:“^(f) — l)dtd2/2. (3-11) 


ry2 

N{y2)= -l)dt 

Jo 


0 JO 
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3.1 Infinite-order homogenised equation: asymptotic approach 

Denote by Hj = Hj(x, y) the coefficient in front of j = 0,1, 2,in the expansion (|3.4I) : 

Hj := curlx{^curlx(./D^-^^V^“^curlxv)| 

+ curly{ j 4 curlx(^^-^'''^^V4curlxv) | + curlx^j curly V^curlxv)} 

+ curly {A curly (.yK*''^’''^^V4''“^curlxv)}, j = 0,1,2,.... (3.12) 

We aim to write Hj(x,y) in the form /i('^'’'^)(y)vi^^curlxv(x), by commuting all x-derivatives through 
to the right-hand sides of the tensors in expressions (j3.12p . To this end, we introduce tensors 

and of order j + 1, such that the following operator identities hold: 

./#(^)(y)Vx = curlx=.T'(^"^)(y), j = 2,3,..., (3.13) 

.i^'^'^^(y)Vx = curlx(^^(y){curly^(^“^)(y)-F./#(^"^)(y)}^, j = 2,3,... . (3.14) 

The expression for Hj, j = 0,1, 2,... is now rewritten as 

H,(x,y) = fi(^'+2)(y)Vi+icurlxv(x), 

fi*''^“'"^^(y) := curl|^(y) curl^^-^“'"^^(y)| -|-curl{74(y)./#*''^“'"^^(y)} -|-jSf^-^^^^(y). 
Summarising, the left-hand side of (12. ip takes the form 

OO 

curl{i^(x)curlu^(x)} = ^ e^/i^^+^)(y)V4+^curlxv(x)|y^^^^. 

1=0 

By analogy with the matrix we require the tensors j = 0,1, 2,... to be independent of the fast 

variable y. The resulting system of recurrence relations 

curl{74curl.yl'*''^“'"^^| = —curl{j4.^^-^^^^| — _|_/i(i+2)^ j — q, 1,..., (3.15) 

allows one to determine uniquely the tensors j = 0,1,..., by virtue of the following 

statement. 

Lemma 3.1. Let be a tensor field whose components are differentiable and Q-periodic. Furthermore, 
assume that A is a positive definite Q-periodic matrix field. Then in order that the equation 

curljjdcurl.yL*''^^} = (3.16) 

have a divergence-free, Q-periodic solution whose elements have zero mean over Q, it is necessary 

and sufficient that 

divF(^) = 0, (f(^)) = 0. (3.17) 

Proof. The necessity follows by taking the divergence and average of both sides of equation (j3.16p and 
noting that both A and are Q-periodic. Conversely, due to the fact that the differential expression 

in the left-hand side of (I3.16h defines a self-adjoint operator in the space of divergence-free tensor fields 
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with zero mean, it suffices to note that the right-hand side of (I3.16P is -L^((5)-orthogonal to all elements of 
the kernel of the left-hand side, namelj{3 

• Vy?) = 0, = (3.18) 

for all smooth order j — 1 tensor fields (p and order j constant tensors The proof is concluded by noting 
that (I3.18P follows immediately from (|3.17l) . □ 

By Lemma [3Tl the equation (|3.15p is solvable if and only if = 0, 

or equivalently It follows that 

/i(-^^^)V'^“i"^curl v(x,e) = v(x, e) = curlcurl -|-./#*^'^'''^)})V-^curl v(x, e)^ . 

The “infinite-order homogenised equation” takes the form 

OO 

curljh^^^curl v(x,e)} -|- g-^cur l{h*^'^’''^^V-^curlv(x,e)} = f(x), x G T, (3.19) 

i=i 

where 

/j(i+2) _ ^A|curl^(^+^) j = 0,1,2,..., (3.20) 

are the tensors of “higher-order homogenised coefficients”. They are related to the tensors via the 

operator identity curlx/i^-^"’"^^ = In the index notation, the formula (I3.20h reads 

Remark 3. For the example ^(y) = a{y 2 )I, we have 

/(a-i)-^ 0 0 \ 

/jW = /^jcurloT'^^^ + -^}') “1 ^ I • 

^ ^ V 0 0 (a-i)-V 

Remark 4. Using index notation, the relations (j3.13p - (|3.14p read 

i = (3.21) 

j = 1,2,..., (3.22) 

where €ijk is equal to 1 if ijk G {123,231,312}, to —1 if ijk G {132,213,321}, and to zero otherwise. 

Substituting (|3.3I) into the homogenised equation (I3.19P , we find that the coefficient functions satisfy 
the following sequence of recurrence relations: 

curljh^^^curlvo} = f, curljh^-^^^^V-^curlvfc} = 0, f = l,2,..., (3.23) 

j+k=l, );,fcGNU{0} 

divvz = 0, / = 0,1,.... (3.24) 

1 jf pU) qU) ^j.g tensors of order j, then ■ G*-^! := 
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Proposition 3.1. The matrix is symmetric and positive definite. Further, for any given right-hand 
side f € [C'j^j,(T)] nAf(T), there exists a unique solution sequence v;, I = 0,1,2,..., for 113.23\} - li3.24\ ) such 
that all elements of the sequence have zero mean over T. 

Proof. Using (13. 8j) and the formula (I3.9jl for the homogenised matrix we write 

/i(2) _ ^^(curl ^'■^^+/) + (Acurl +^)curl = ^74(curl (curl ^ (3.25) 

In view of the fact that the matrix A is symmetric and positive definite, the expression in the right-hand 
side of (I3.25h is also symmetric and positive definite. Applying the Lax-Milgram lemma to each equation in 
(I3.23|] . in conjunction with the corresponding equation in (j3.24p . yields the second claim of the lemma. □ 

In Section [3.31 we provide a justification of the above formal argument, by showing that the remainder 
of the truncation of the double series (I3.2p - (l3.3p is close to the solution of the original problem in a certain 
sense. Before doing so, we briefly discuss the recurrence relations that determine the tensors These 

relations arise by considering the divergence of (13.211 . 


3.2 Recurrence relations for = 2,3,... 

Taking the formal divergence of the asymptotic expansion (13.2p and using the condition div = 0 yields 


^e^|(^(^)(y)Vj)v(x,e)) 

i=i 


3,.y.+ -kdiVx(^(^)(y)Vj) ^curlxv(x,e)) | 

oo 

+ + (jr(^-i)(y)Vi-V(x,e))^^^^^} = 0. 

i=i 


Comparing the terms with equal powers of £ yields the system of recurrence relations 

(jr(')V,v) = 0, (3.26) 

, 2 ,... . 

(3.27) 

Note that the Q-periodic solution to p3.26p with zero average over Q is identically zero. Hence, the 

first non-trivial tensor in the sequence is the third-order tensor which satisfies 

(jr(2)v2v)^^^^^ = -div4^(i)curl,v), 

for an arbitrary vector v. Substituting the expansion p3.3p into the system p3.27p . we obtain 


+div4.yK«V(r'curl,v) = 0, / = 1 




E 

j+k=l 

ieN, fceNufo} 


,yiyi 






+ diVx(./K<^’Vi ^curlxVfc) I ^ 0, 


; = 1,2,... . (3.28) 


We use the set of equations ^3.2811 to establish a bound on the divergence of the solution and, subse¬ 
quently, on the remainder of the asymptotic series p3.ip , see Section 13.31 
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3.3 Remainder estimates 

The justification of the above formal procedure is given by the following theorem. 

Theorem 3.1. For all K consider the remainder e) := u^(x) — u('^)(x,e), where 

K 

uW(x,e) := vW(x,e) + ^e^'{Vy(.^(^)(y)V4vW(x,e)) + Vx(jr(^-^)(y)Vi-VW(x,e)) 

i=i 

_l_ ^0') (y) (curlxV^”^^ (x, e 
K 

v(^)(x,e) := ^e^Vfc(x). 
k=0 


y=x/£ 


(3.29) 

(3.30) 


Then the estimates 

(i) ||curlR(^)||^2^.j,^ < 

(ii) ||divR(-^)||^_i^^^ < 

(iii) VM |(rW)^| 

hold, where the constants are independent of e but may depend on the function f. 

Proof, (i) We evaluate the operator on the left-hand side of (|2.in on the function : 

K-2 

curlji^curlu^-^)} = ^ e^curlxjh^^'+^^curlxV^-^)} + e, K) , (3.31) 

j=o 

0i(«;e, ii') := e ^curlx{^curlx{Vy(jT^*''^^V^«)} -|-c/h^'^^V^“^curlx»}^ 


+curly{AcurlxVy(jr(^)V^.)} + {curly (^.^(■^+1)) + curlx* 

Substituting the expression for from (13.301) into (|3.3ip . we obtain 

K-2 K 

curlj^^curlu*-^^} = £'^'*~^curl{h^'^'*~^^V'^curl v^j 

j=0 k=0 


K 


+ e^ ^ ^e^©i(vfc;e,if) 
k=0 


i + e^ ^&2{x,e,K), 


2K-2 K 

& 2 {x,e,K):= e^+*^-'^+icurl{/i(^+2)V'’curlvfc}+ ^e^0i(vfc;e,K), 

j-\-k=K—l k=0 

where we use (j3.23p . It can be shown that 

|02(x,e,if)| xeT, 
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with an e-independent constant > 0. Recalling (|2.1I) . we obtain 

curl{A^(x)curlR*^'^^(x,e)} = —e^“^© 2 (x, e, iT). (3.32) 

Further, taking the scalar product of both sides of equation (|3.32p with integrating over T, and using 

the Cauchy-Schwarz inequality on the right-hand side yields 

j i^curlR(^) •curlR(^) = J ©2 • R^-^^ < e-^"^||©2||[z,2(T)]3||R^^^||j^2(T)]3- 

Since the matrix is positive dehnite, it follows that 

z2||curlR(^)||j^2^^^j3 < e-^"^||©2||[L2(T)]3||R^^^||[i2(T)]3, 

and using the Poincare-type inequality from Appendix A yields 

i/||curlR(^)||j^2(T)]3 < e^-^|T|||©2||[L2(T)]3(||curlRW||j^2(^^j3 + ||divRW||^2(^)). 

It will be shown in the proof of (ii) that ||divR(^^||^ 2 (T) ~ Hence 

||curlR(^)||j^2(T)]3 < |T|||©2||[L2(T)]3Z2"^e^"^||curlR(^)||j^2(T)]3 + 0(6^^“^), 

which implies the claim. 

(ii) Note that divu®" = 0 and hence divR^^^^ = —divu*^^\ Denoting by U/ the coefficient in front of 
in the combined sum (I3.29p ^ (|3.30p : 

K ra\'a{l-j,K} 

Uo:=vo, Vi-=Y1 Y1 {Vy(^(^)vivfc)+Vx(^(^-^)vrVfc)+^(^)vrWlxVfc}, 1 = 1,2,..., 

j=l k=0 

(3.33) 

we write 

2K 

-divR(-^)(x,e) = ^{e'diVxUz + e^“MivyUj 
i=i 

Taking into account the recurrence relations (I3.28p . we obtain 

2K 2K 

-divR(^)(x,e) = e^03(x, e, AT), 03(x, e, AT) := < ^ e^“^diVxU/^ e^“^“MivyUz 

l=K l=K+l 

Note hrst, see (I3.33p . that U; is a finite sum of terms of the form C/(y)P(x), for some tensors U{y) with 
elements in H^{Q) and R(x) with elements in C'“j,(T). Further, since 

= ediv 

the two-scale form of 03 is a finite sum of terms of the form U (y)P (x), for some tensors U (y) with elements 
in A^(Q) and R(x) with elements in C'°°(T). Finally, we use the following statement, which is a version of 
the theorem in m Appendix C]. 


1 y=x/£ 
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Lemma 3.2. Let M(x/e) he a periodic tensor of order j whose components have zero average on Q and 
let g{x.) be a smooth, periodic tensor of order j — 1. Then there exist positive constants r G N, such 
that 


M(x/e)^(x)dx 


< 


Vr G N. 


(3.34) 


It follows from the above lemma that the L^(T)-norm of ©3 is bounded by an e-independent constant 
6 * 2 ^ \ and hence 

||divR(^)||^2(.jj,) < e'^l|03||i2(T) < 

(hi) Integrating (I3.29P over T and using Q-periodicity of and T-periodicity of v yields 
f dx = f o4^*'-^^(x/e)V-^“^curl v^^^(x) dx. 

Jt ^ Jt 

The result follows by applying once again Lemma 13.21 
Corollary 3.1. The estimates 

(t) < cf) > 0, 

curl' ' 

hold, where and are e-independent but may depend on the function f. 

Proof. The proof is an immediate consequence of the theorem and the inequality (j6.25p . □ 

Remark 5. The above method of truncating (I3.19P may lead to non-elliptic higher-order problems. An 
alternative approach, which is free of this limitation, is discussed in Section [H 


Ir(^)II ^ r'^K) K-i 


□ 


(3.35) 


3.4 Example: laminate with two layers per period 

Suppose that A(y) is given by 

A{y)=a{y2)I, a{y2) = \ / <1 

P 02, T < 2/2 S 1, 

where 0 < /i < 1, and oi, 02 are positive constants. In what follows, we determine the leading order terms 
in the asymptotic expansion (|3.19l) 

curl{A^curl u^} = curlxjh'-^^curlxv} -|-ecurlxjh^^^VxCurlxv} -|-e^curlx{/i^^^VxCurlxv} -|-0(e^), 

for the case above. For the term of order 0(1), the related tensors are given by 

/O 0 -N\ 0 0 \ 

= 0 0 0 , = 0 ° ’ 

\iV 0 0 / \ 0 0 (a-i)-7 

where N = N{y 2 ) satisfies the differential equation —{aN')' = a' see Remark [2l and 

( 0 “^)“^ = [cxf^h + o^^(l - /i)}~\ (o) = (oih -F 02(1 - h)). 
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For the term of order 0(e), we have 


< = < 


M, ijk = {123}, 
—M, ijk = {321}, 
L, ijk = {132}, 
-L, ijA: = {312}, 
0 , otherwise. 


ijk = {112, 332}, 
ijk = {211 
0 , otherwise. 


hfk = { -b, ijk = {211, 233}, 


where 


{aM'Y = (aN)', {aL')' = (a) — a, a = —{aL'), b = (aN). 

Finally, the relevant tensors for the term of order O(e^) are given by 


(3.36) 


^hjkl 


= 


r p, 

ijkl 

= {1232}, 






-p, 

ijkl 

= {3212}, 


c. 

ijkl 

= {1212, 

3232}, 

Q, 

ijkl 

= {1223}, 


d. 

ijkl 

= {2121, 

2323}, 

-Q, 

ijkl 

= {3221}, 

'hjkl ~ 

e. 

ijkl 

= {2112, 

2332}, 

R, 

ijkl 

= {1311,1333}, 

/, 

ijkl 

= {1111, 

1133,3311,3333} 

-R, 

ijkl 

= {3111,3133}, 


1 0, 

otherwise. 


[ 0, 

otherwise. 






-{aP'] 

' = 0 

-|- oiL^ -\- {^(xLy ^ 

-{aQ')' = 

{aM)', 

-{aR')' = 

= b-aN, 


where 

— {a 

c = —{aP' + aL), d={aQ'), e = {aL), f = —{aR!). 

It is shown by direct calculation that 

a = b = 0, c= —d = — 12 “^/^/|(q:^^/i + a^^h) ^(/3i — l )(/32 — 1 ), 

e = + aY^l2)iai - 02)^, / = 12 ~^ ^^2)~^(1 - /3i)(l - /32)(a^^/i + a)“^^2), 


where /i + ^2 = 1, and /3i = 02 ^ai = (32^ ■ The homogenised equation takes the form 


3-1 


/(q:“^)“^(u 2,12 - 1^1,22) - (a)(^^l,33 - 'y3,13)\ 
(q:“^)“^(u3,23 - 1^2,33 “ ^^2,11 + ^’1,12) 
\(a)(?^l,13 - ^’3,ll) - (a“^)"n^3,22 - '^ 2 , 23 )/ 


+ £^ 


^ /(ns, 1223 — ^1,2233) — e(uiqi 33 — usqus + 111^3333 — 1134333) 

0 

\—/(ns,1122 — ni,i223) + e(ni,iii3 — ns^m + U14333 — 1134133)^ 


+ 0{£^) = f 


Remark 6 . It was shown in | 12 j that in the case of a scalar equation, all terms with odd powers of e are 
absent from the corresponding infinite-order homogenised equation. The above two-layered case provides 
an example in which there is a non-trivial term of order O(e^) in the homogenisation procedure for the 
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Maxwell system. It is calculated that 


= 


A^i, ijklm = {12232}, 

—Ni, ijklm = {32212}, 

N 2 , ijklm = {13121,13323}, 
-N 2 , ijklm = {31121, 31323}, 
N 3 , = {13332,13112}, 

-Ns, ijklm = {31332, 31112}, 
Ni, ijklm = {12311,12333}, 
-Ni, ijklm = {32111, 32133}, 
Ns, ijklm = {12223}, 

—Ns, ijklm = {32221}, 

Nq, ijklm = {23212}, 

-Nq, ijklm = {21232}, 

Nt, = {23111,23133}, 

-Nt, ijklm = {21311, 21333}, 
0, otherwise. 


= 


hi, ijklm = {12212,32232}, 
f 11121, 11323} 


h 2 , ijklm = 
hs, ijklm = 


{33121, 33323] 

| 11112 , II332I 

{33112, 33332] 
hi, ijklm = {21212,23232}, 
hs, ijklm = {21221,2322^}, 


he, ijklm = 


h-j, ijklm = 

0, otherwise. 


f 21111, 211331 
{23311, 23333] 

h2iii, 12133I 

{32311, 32333] 


so that the term in question is given by 

{hi — h2 + hs, — /l7)(ui,22233 “ '^3,12223) + (^3 — hi + /l6)(uiqi233 “ '^3,11123 + f1,23333 “ f3,12333) 

0 

y—{hl — /i 2 + /15 — /l7)(t;iq2223 — 'f3,11222) “ (^3 — hi + hQ){vi^i2333 “ ^3,11233 + 1^1,11123 “ '^3,11112)-^ 

The functions Ni, i = 1,... ,7, are functions of 1/2 only, which satisfy 

-{aN[y = {aPy + aP'+ aL + c, -{aNyy = aM + d, -{aNf^)'= aL - e, 

-{aN'i)'= {aPy+ aR'+ f, -{aN'^)'= {aQ)', -{aN'o)'= aP'+ aL + c, -{aN'j)' = aP + f, 

and the constants p, i = 1,... ,7, are given by 

hi = —{aN[ + aP), /12 = —{aN^), /13 = —{aN^}, 

hi = (aP), hs, = {aQ), h^, = {aP), hj = —{aN^ + aP). 

It is shown directly that hi = —hs,, h 2 = —hj and /i 3 = —hg, hence the above order 0{e^) term takes the 

^—hi{vi^ll233 — 1^3,11123 + 1^1,23333 “ 'f3,12333)^ 

0 

hi{vi^l2333 — 1^3,11233 + 1^1,11123 “ 1^3,11112) / 

4 Variational approach 

Consider the minimisation problem 


form 


mm 


-A^ curl u • curl u — f • u) =: /(e, f), 


(4.1) 


where the minimum is taken over T-periodic functions with zero average and zero divergence. Clearly (EH) 
is the Euler-Lagrange equation for the functional (14.11) . Therefore, one has 


IieA) = --. 


curl u^ • curl u^ = — 


1 


f u^ 


(4.2) 
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4.1 Variational asymptotics 

Similarly to (SH), the solution to the homogenised equation (13.9p is the solution to the minimisation 
problem 

min J ^-h*'^^curl V • curl V — f • =: /o(f)) (4-3) 

over all divergence free, periodic functions on T. It is well-known, see e.g. [9], that the functional m 
converges to the homogenised functional (j4.3p as e —)• 0, i.e. for any function f one has /(e,f) —^ lo(f) as 
e —)■ 0. The following result generalises this fact to all finite orders in e. 

Proposition 4.1. For any function f € [d^rC^)]^ positive integer K, there exists a constant Ck 

such that 

I(e,f) + i [ f(x).v(^)(x,£) dx 
^ Jt 

Proof. Using Theorem 13.11 we write 

= [ f.uW+ [ R(^)(x,e)dx 

^ JT JT 


< Cks^. 


- [ f(x) • vW(x,e)dx+ V [ ^(^■)(x/e)Fjfc(x)dx + [ R(^)(x,e)dx, (4.4) 

2 Jt jt Jt 


^<j<K 

0<k<K 


where 


R^'^^(x, e)dx 


< Ck£^, Ck > 0 , 


and Fjk := Vx ^curlxVfc * f are infinitely smooth and T-periodic. As the tensors are periodic and 

have zero mean, the integrals 




[ .yU(^)(x/e)T;fe(x)dx, 
JT 


1 < J < at, 0 < k < K, 


go to zero as e —> 0 faster than any power of e, see Lemma 13.21 This observation and (|4.4I) immediately 
imply the claim. □ 


The above proposition can be interpreted in the sense that 


OO ^ « 

I(e,f) Ve"4(f), 4(f) / f-vfc A: = 0,1, 2,... . 

k=i ^ 


(4.5) 


4.2 Infinite-order homogenised solution 

In this section, the effect of the rapid oscillations in the tensors in the asymptotic expansion (13.2p 

is removed by following the translation averaging strategy of [12]. More precisely, for all G Q, denote 
A<^(y) := A(y + C); y € Q, and consider the equation 

curl{ A^(x) curl u(x)} = f (x), x € T, (1-6) 
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under the same assumption for the right-hand side f and subject to the same conditions on the solution u 
as in the case of (j2.ip . For all C) the equation ()4.6I) admits a unique solution u^. Consider the average of 
with respect to i.e. 


u^(x) := (u^(x))^ ^ ^ 


(4.7) 


IQ 


This averaging with respect to C is analogous to “ensemble averaging” for the family ()4.6p when the 
underlying probability measure is uniform over Q. The next proposition shows that translation averaging 
eliminates the effects of oscillations due to the tensors and results in the validity of the asymptotics 

for the average u^. 


Proposition 4.2. For a given function f, the series i3.3\} provides asymptotics for the function in the 
sense that for any positive integer K, there exists a positive constant Ck such that 


<CKe^. 


Proof. Note that for all Q Q, an expansion similar to (I3.2p - (l3.3p holds for the solution of (I4.6p . where 

and are replaced by the “translated” tensors := and '■= PF^ii\- + C), 

respectively. In particular, by analogy with the first estimate in (j3.35D . one has 


I £ iK)\\ 

Pc p ll[L2(T)]3 


< Cke^, 


Ck > 0 , 


where is a truncation of the asymptotic series for similar to (I3.29|] . with all y-dependent objects 
replaced by their (^-translated versions: 


uf\^,e) = vP)(x,e) + J]p{Vy(jr/^)(y)VivP)(x,e)) 


i=i 


+ Vx ((y) vr VP) (x, e)) + (y) vr 1 (curl^vP) (x, e 


U), 


(4.8) 


Next, we note that 


uPx)-vP)(x,e)= [ (u-^(x)-uP)(x,e))dC, 

JQ 

since the integrals of the y-dependent tensors under summation in (14.Sp clearly vanish. Therefore, one has 


^ |uPx) - vP)(x,e)|Vx < |u"cV) - uP)(x,e)|dCydx 




/ .,2K 


C'k>0, 


Q JT 


as required. 


□ 


4.3 Higher-order variational problems 

By analogy with (j4.1l) . we define the “translated” energy functional: 

/^(e,f) := niin 

u(x) 


curl u • curl u — f • u 


(4.9) 
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and the “averaged” functional 


J(e,f) := [ I^(e,f) dC- 

JQ 

Considering (^-dependent trial fields in (14.91) and changing the order of (^-integration and minimisation 
yields 

/(e,f) = min .Ee(u, f), (4.10) 

u(x,C) 

where 

.Ee(u,f);= J j ^-^^(x)curlu(x) • curlu(x) — f(x) • u(x)^dxd^. (4-11) 

Clearly, the variational problem (I4.10p has as its minimiser the function u(x, (^) = u^(x), where u^(x) is 
the solution of (14.6p . Further, recall that is represented by the series 


OO 

u-^(x) = v(x) + ^e^{Vy(jr/^')(y)Viv(x)) + Vx(jr/^“'^(y)Vr'v(x)) +^/^)(y)Vr'curlxv(x)} 


i=i 


y=x/e 


For each positive integer K, consider the subset of the set of trial fields in (|4.10p obtained by truncating 
the above expansion: 


Uk-= 


K 

u(x,C) : u(x,C) = v(x) + j;e^{Vy(jr/^')(y)Viv(x)) 


i=i 


+ Vx(=^*''^ ^^(y)V^ ^v(x))-|-^curlxv(x)| for some v|. (4.12) 


Here v runs over the set [C^j,(T)]^ n d:’(T) of T-periodic, divergence-free, smooth vector fields with zero 
mean. Consider a minimisation problem for the same functional ()4.10p over the set Uk- Substituting test 
functions from (I4.12p into the equation (I4.10p yields 


.KK . 

£i+fc)i(hfc)VJcurl vV^curl V - f • v L 

j=0 k=0 ^ ^ 

where is a tensor of order j -|- A: -|- 2 defined by 

~^U,k) ^^(curl^(^'+i) +.^(^’+i))(curl^(^+^) 

or in the index notation: 

Cl,™ = ('4.4curl..r(J+‘> 

Whenever u and v are related via the expression in the definition of Uk, we set 

EK{^,i,e) := Ee{u,i), 


(4.13) 


(4.14) 


so that 

min£'i^(v, f, e) = min .^^(u,!), (4.15) 

V u(x,C)€!7k 

where v € [C^j-CIT)]^ C fF(T). Clearly for all e, f, the functional Fiii'(-,f,e) is convex, as a result of C- 
averaging of the convex functional The next proposition, which is similar to [12p Proposition 3], 

shows that the minimiser vk in ()4.15p is the best variational choice of a truncated approximation. 
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Proposition 4.3. For K >2 and all functions v one has 

EK{v,i,£) > EKivK,f,£) > 

Moreover, the estimate 


-Ex(vx,f,e) - < cke 


2K 


(4.16) 

(4.17) 


holds for some eonstant ck > 0. 


Proof. The inequalities (|4.16l) follow from the minimising property of vk and the fact that Iii:(e,f) > 
/(e, f). To obtain the estimate ()4.17p . we substitute u^(x) = u^^^(x, e) + R^^^(x, f, e) into (14.ip and 

integrate by parts. The remainder satisfies bounds analogous to those proved for R*^^) in Theorem 

EH hence 

0 < ^ ckE^^. 

Integrating the last inequality with respect to ^ G Q and using the minimising property of v again, yields 

(HTTP . □ 

Note the following alternative formula for _ Uging integration by parts, for j = 1, 2,..., we get 

(7l(curl + .^(^■+i))curl f) = G [C'^,(g)]^ {(f) = 0, 


hence 

where 


h^hk) ^ (^(curl 

Ar(k+P\ _ (^{j+P yik+i) 

'' -^sn.-.h+i ^ 


(4.18) 


'll...lj+k+2 Sll...lj+i Slj+2---lj+k+2 

4.4 Infinite-order variational homogenised equation 

The weak form of the Euler-Lagrange equation for the problem (I4.15P reads 

curl (^V*^curlvii-+V-^curl Vi^-V^curl—f-(^1 =0 V(^G [(^^^.(T)]^. (4.19) 
j,k=o ^ 

For K = 0 this gives the identity 

'-^(°’°)curl(^ • curl Vo + -^(°’°)curl vq • curl - f • = 0 V if & [C“,(T)]^ 


or in the differential form 

curl|/i*'°’'^^curl vq} = f, 

where as follows from (14.1411 . Proceeding in a similar manner, the Euler-Lagrange equation 

for the case K = 1 reads 

curljL^'^’^^curl vi} -|- ecurlj^^Vcurl vi} -|- curljL^^’^^Vcurl vi} = f. 


71 ._ 1/7(04) _ r(0,l) _ r(l,0) r(l,0) 

^ijk 


^ijk ^kji ^jik ' jki 


}■ 


(4.20) 
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Consider the case when formally X = oo in (14.191) and denote by the corresponding minimiser. Notice 
that 

min .Eoo(v, f, e) = min f) = I(e, f), 

V u(x,C) 

where the minimiser u(x, (^) = u^(x), and therefore = u^. After a series of manipulations, the related 
Euler-Lagrange equation is written as 


curl{h(°’°)curlv(°°)} + ^e”curl{h"V"curlv(°°)} = f, 


n=l 


where (c/. (14.141) 1 


/i” ■ 


-5 E {(-1)^'^ 


ij+li2-"'ijil'ij+2---ij+k+2 


+ i-irh. 


kZ(j,k) 




(4.21) 


(4.22) 


j-\-k=n, 

j,/cGNU{0} 


In Section O we show that the “asymptotic” infinite-order homogenised equation (|3.19p coincides with the 
“variational” infinite-order homogenised equation (|4.21l) . 

The next result shows that the minimiser vk of Eii-(v, f,e) is an approximation of order 0(e^^) to 
the infinite-order homogenised solution. 

Proposition 4.4. Let K >2 be a positive integer. Then for any vector function f, there exists a constant 
ck > 0 such that 

II""-v^II[l2(t)]3 <CKe^- 

Proof. The proof follows the argument of m Appendix D]. Denote by G Uk the vector associated 
with yk by (14.121) . and let R^(x) := u^(x) — u^(x,e). Notice hrst that 

1 


£’K(v/f,f,e) = E£(u^,f) = f E£_^(u^,f)dC = [ [ ( 

Jq Jq Jt \ 


-A^ curl • curl u^ — - curl u^ • curl 


- curl u^ ■ curl u^ — f • u^ ) dxd(^ 



Q JT 


-A^ curl R^ • curl u^ -|—A^ curl R^ • curl R^ — f • u^ -|- f • R)^ | dxd^ 


Since solves (14.6h . integrating by parts in the last expression we obtain 

Ei^(vi^,f, e) =/(e,f) f / ^Af curl R^ • curl R^dxdC- 

Jq Jt J 

The fact that the matrix A^ is positive-definite implies 


A5 curl rL • curl R^dxdC > v 


Q JT 


QJT 


|curlR^|^dxdC 


> n 


TJQ 


[curlR^l^d^dx > f |curl (u^ — vx)|^. 

JT 


(4.23) 

(4.24) 

(4.25) 


Combining (I4.23p . p4.25p . and p4.17p yields 

|curl (u*^ — 'vk)\‘^ < 2cj^zj“^e^^. 


Finally, noticing that divu^ = divvj^ = 0 and using (|6.25p . we obtain the desired result. 


□ 
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5 Tensor analysis of the infinite-order homogenised equations 


Here we prove that the infinite-order homogenised equations (I3.19P and (14.211) are equivalent. To this end, 
we introduce a symmetrisation operation, as follows. 

Definition 5.1. For n > 1, the symmetrisation of a tensor of order n -|- 2 is a tensor of 

the same order defined by 

1 

T'iki...knji 

{ki,k2,...,k„) 

where the summation is over all permutations of the indices ki,k 2 , ■ ■ ■ ,kri 


hi(ki...k„)j ■ ^ hi 


Notice that in the case of third-order tensors, the above operation leaves the tensor unchanged. In 
order to account for this special case, we show separately that the third-order tensors and coincide. 
To this end, notice first that 

= + /)^^(curl.yF(2) + . 

Further, making use of the definitions (I3.2ip and (|3.22p and relation (I4.18P yields 

h% = ^ ^( 0 , 1 ) _ ^( 1 . 0 )_ 

Hence, by virtue of p3.2np and (I4.20p . one has 


_ 1(0,1) _ 1(0,1) _ 

“ "‘ijk '‘■kji 




7(0,1) _ 7(0,1) _ 7(0,1) r (0,1)1 
"'ijk "kji "kji "ijk j 


= 0 . 


Remark 7. It can be shown that for a third-order tensor h = (hijk) the identity curI{/iVcurI v} = 0 
holds for all v if and only if 

oi, fjA: G {122,133}, 

02 , fjA; G {221,331}, 

03 , ij/c G {211,233}, 

04 , ij/c G {112,332}, 

05 , fjA: G {311,322}, 

06, fjA: G {113,223}, 

01-1-02, ijA; G {111}, 

03 -|- 04 , ijk G {222}, 

05 -|- 06 , ijk G {333}, 

0, otherwise, 

for some constants o;, I = 1, ..., 6. It follows from the above that in the case of the tensor h = — h^, 

the constants o; vanish for all 1. 


hijk — ^ 
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Returning to the case of an arbitrary order, consider the expression for the tensors appearing in the 
variational approach (c/. (j4.22l) l: 

h"-, , ^ I .1 

2 / ^ ' kph\...hp—±ikp-^\...hp-^-qj ^ ^ hph\...hp—\jhp-^\...kp-^qij 

p+q=n, 
p,qeNU{0} 

Symmetrising the above expression with respect to ki, /c 2 , kn gives 


(5.1) 




277,1 / ' kpkx • •-kp—xikp-^i^.. .kp-^qj V f kpk\...kp—\jkp^\...kp^q'i 


(fcl,...,fcn) P+q=n, 

p,ijeNU{0} 


}• (5.2) 


The symmetry property 


h) 


(p^q) 




ih\. . .fcp/Cp-j_l . . .fcp_|_gjl /Cp_j_l . . .fcp_|_gjl 2 /ci . ..kp 


implies 


- V V 

^ ' /Cpfcl...fcp_l 2 /Cp-j_l...fcp+gJ' 


(t E E 


kpki.. .kp-ijkp-^-x.. .kpjf^qi ^ 


(ki,...,k„) p+g=n, 
p,qeNU{0} 


(ki,...,kn) P+q=n, 

p,qeNU{ 0 } 


and hence 


In _ J_ r_T'\Ph(P’9) 

'‘■z{ki...kn)j 


E E (-1)"^. 


kpk \ . ..kp—\ikp^ \ ...2p-|_gj' 


(5.3) 


(ki,...,kn) P+q=n, 

P,<jeNu{o} 


As is shown above, the third-order tensor is expressed as a combination of the tensors 

Similarly, the tensor is a combination of the tensors ..., as follows: 


? (n+2) _ lAnfi) _ r (n-1,1) 7 (n-2,2) _ r (n,-3,3) , 

^iki...knj ^ki...knji ‘'k2...knjkii ' ^k3...knjk2kii ''ki...knjk3k2kii ' • • • 




(0,n) 

jknkn—1 ...k\i 


The related argument makes use of the equations (|3.21l) - (|3.22l) and the relation (|4.18p . Symmetrising the 
above expression and using symmetry properties of the tensors , we obtain 


h 


(n+2) 


i{ki...k„)j V ^ki...knji ’^k2...knjkii 

(fcl ,...,kn) 


Anfi) , |'_iNn7(0,n) 1 

'll., u i, -I- -I- l, 1) ll-jk^k„-i...kii} 


1 


n! 


Z { 

(fcl ,...,kn) 

Yhki...k„j ^ '‘'kAn-l-klij j 


{k\^...^kn) 


2; E E (-i)”'" 


(p-q) 


kpk \.. ./Cp_ 1 .. .ip^qj 2(^1 ■ ■ -kAj 




(5.4) 

(5.5) 

(5.6) 


{ki,...,k„) P+q=n, 
p,q€NU{0} 


Hence, the following result is established (c/. [121 Section 3.3.2] for the scalar case): 

Theorem 5.1. For the tensors and h” given by lld.2(j\) and {4-22^ respectively, one has 


f (n+2) _ Tn 

i(ki...kn)j h(fci...fcn)j' 


In particular, since each of the infinite-order homogenised equations i2.iy\) and ^-21^ is unchanged under 
this symmetrisation process, these equations coincide. 
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6 Homogenised constitutive laws for the system of Maxwell equations 


6.1 Maxwell equations and qnasistatic approximation 

Written in time-harmonic form, the system of Maxwell equations for the electric and magnetic 
fields and electric displacement and magnetic induction in an e-periodic medium is 


curlE"^’^ = -iwB'^’^ 
divD‘^’^ = 0, 
B“'^ = 


curlH“’^ = iwD‘^’^ + 
divB‘^’'' = 0, 


( 6 . 1 ) 

( 6 . 2 ) 

(6.3) 


where w > 0 is the frequency, and we assume that —divJ‘^ = 0. The coefficients /r, e in (|6.3p are the 
matrices of magnetic permeability and electric permittivity, which are subject to the same assumptions as 
the function A at the beginning of Section [2j Denoting generically by A the fields entering (|6.ip ~ ()6.3p . we 
substitute 

OO 

A-’^(x) = ^(iw)^ A^^(x), A^^(x) € M Vx G T, j = 0,1,2,..., (6.4) 

j=0 

and compare the terms with equal powers of the frequency uj. This procedure, which we refer to as the 
“quasistatic approximation”, results in the leading-order equations 


curl{(/i^) ^ curlEf} = — Jo, curl{(e^) ^ curlHf} = curl{(e^) ^Ji}, 
div(e^Ef) = -div Ja, div(/i^Hf) = 0. 


The focus of Sections [2] and [His the first equation in the system (|6.5I) . where we set Ej =: with — Jq =: f. 

Note also that under the assumption that the electric permittivity is given by the identity matrix, one has 
div = 0. 


Remark 8. The form of the equation for the magnetic field is somewhat different to that for the electric 
field. Indeed, the right-hand side of it depends on the permittivity matrix and hence depends on e. A 
modified approach that deals with this feature is described in Section 16.41 


6.2 Energy considerations 

The total energy of the electromagnetic field is the sum of two parts: the electrostatic energy and the 
magnetic energy. For the system of equations (I6.ip ^ (|6.3l) . the corresponding expressions are given by 


1 




■= 

‘^mag • 


1 




( 6 . 6 ) 


where the bar stands for complex conjugation. Using equations (|6.3p - (|6.ip in (|6.6I) . we obtain 

1 


^elec ■ 


2w2 


[ (e^)-^(curlH‘^’= - • (curl - J^), 

Jt 

[ 

Jt 


u‘^’^ ■= 

“mag • 2a;2 


curl E“^’^ • curlE*^ 


(6.7) 


( 6 . 8 ) 


Note that in the modified expression the electric energy depends on the magnetic field and that the magnetic 
energy depends on the electric field. From here on, considerations will be restricted to real vector fields 
and hence the complex conjugation notation will be dropped. 
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Substituting a formal series in powers of oj for D^’^, B'^’^ into (I6.6D - (I6.8D and comparing 

the coefficients in front of (iw)-^, j = 0,1, 2, yields: 


{jf) ^curlEg • curl Eg = 0, 


(/l^)-^curlE^ • curlE^ = / 




H-. 


(6.9) 


Note that due to the fact that curl Eg = 0, the first equation in (|6.9I1 holds automatically. 

In the next section we use the approach developed in Section [2] and Section H] to derive higher-order 
constitutive relations between the leading-order magnetic Hg and induction Bg fields. The asymptotic and 
variational approaches lead to two expressions for such a higher-order constitutive law, which are shown 
to coincide, by a symmetrisation procedure in Section [5l 

6.3 Infinite-order constitntive relations 

In this section, we derive an expression for the leading-order term Hg, given that the expansions for the 
magnetic field Bg and the electric held Ef are known. In the quasistatic approximation we obtain 

curl{(/i'^)“^curlE^} = f, Ef G -T(T), 

where f := — Jq G [C^j-Cir)]^ fl d:’(T). These equations were analysed in Section [2l in particular, Ef is 
written as a double series (|3.2p ~ (|3.3p and the inhnite-order homogenised equation ()3.19p is satished. 

We establish a higher-order constitutive relation between the magnetic held and magnetic induction 
by two approaches. First, we examine the family of problems considered in Section [4.21 In particular, we 
introduce a parameter Q acting as a shift in the microscopic variable, which is rehected in the notation 

E^’^. Notice hrst that since is the identity matrix, the asymptotic expansion for the electric 


-c 


held displacement is simply the asymptotic expansion for E^’^. Hence, averaging over Q with respect 
to ^ yields Df = Ef = v, where the bar notation represents averaging over Q with respect to Q. Further, 
note that 

- /r.ny'l -\T ^0 16) 


and therefore 


— (Bq = (curl Ef’^)^ = curlv, 

OO 

: = = -((/i^)-^curlEf’^)^ = - ^ e^'/i^^'+^iv^curl v 


3=0 


by 1I4.8I1 


J]V^B^. 


i=o 


( 6 . 11 ) 


i=o 


The resulting constitutive law is of the form Hg = (/r®®)“^Bg, where (/i®®)“^ is the inverse of the “effective 
permeability operator”. 

In [12], a higher-order stress-strain relation is derived via a variational argument by considering the 
corresponding elastic energy functional. Applying the same approach to the magnetic energy functional in 
(lOD . an alternative expression for the higher-order constitutive law is derived. Results analogous to those 
of Section m hold, in particular: 




curlEf’** • curlEf’*’ = — 


£.C 


1 


■‘1 


1 


1 


« OO 

/ f-Ef’^ ^^£"^4(1), 

k=0 


4(f) = -x / f ■ Vfe, 


A: = 0,1,2,... 
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( 6 . 12 ) 


7(e,f) 


1 

2 



V^curl vV^curl 


V, 


where is the tensor given by (|4.22n . Integrating by parts in (|6.12p and using (|6.10l] yields, after an 
appropriate rearrangement of indices: 


vf 

^mag 



oo 

n=0 j-\-k=n 

j,A:GNU{0} 


(6.13) 


where 

hli2...b+fe+2 ■ ^ii+li2---ijij+k+2ij+2---ij+k+lh' 

Using the symmetrisation procedure described in Section [5] it is confirmed that the expressions (16.110 and 
(16.1311 for the magnetic field coincide to all orders. 


6.4 The magnetic field equation 

In this section we assume that the permeability is the identity matrix and that the permittivity is 
Q-periodic, symmetric and uniformly elliptic. We consider the quasistatic approximation of Section (|6.2p . 
For a given vector function Ji € [Cj^^C^)] ^(T), we look for a solution to the problem 

curl{(e'^)“^curlHf} = curl{(e^)“^Ji}, Hf € A:’(T). (6-14) 

The theory discussed in Section [2] applies, subject to a modification to the asymptotic expansion (13.2p : 


(x) = w(x, e) + ^ { Vy (.^ 7) (y) V7w(x, e)) + (y) V(x, e)) 

+ =^(''^(y)Vi-i(curlxw(x,e) - Ji(: 


i=i 


X 


y=x/e 


, (6.15) 


where w G [C“j,(T)]^ is a divergence-free vector field written as a 


series 


w(x, e) = ^ e*^Wfc(x), x G T, (6.16) 

fc=o 

and are tensor fields of order j + 1, j = 0,1,... , = 0. The equations of Section [2] 

hold with replaced by replaced by and 7^7 replaced by where ^7) = 

^(i+i)Vx = curlx.;T(-^7 In particular, the homogenised equation of infinite-order has the form 

OO OO 

^^e-^curl{7^-^^^^V'^curlw(x,e)| = e'^curl|fc^-^'*~^^V-^ Ji(x)}, (6-17) 

i=o i=o 

^(j+2) _ (r^{curl.^(^+^) +.^(^'+7}). 

A version of Theorem l3.1l on iustihcation of the asymptotic procedure is proved, with similar error estimates. 
When considering the variational approach developed in Section [H the argument is modified as follows. 
Consider the expression for the electric energy given by the first formula in p6.6l) . Applying the ideas 


23 













discussed in Section 16.21 we first apply the quasistatic approximation to the energy functional and then 
write it in terms of the magnetic field which yields 

/e^-D^= /(e*')-!(curlHf-Ji) • (curlHf-Ji). (6.18) 

JT Jt 

It is straightforward to see that (I6.14p is the Euler-Lagrange equation for the minimisation problem 

mini f (e^)“^(curlH — Ji) • (curlH — Ji) =:/(e, Ji), (6.19) 

H 2 Jt 


Analogous conclusions to those given in Section 0] hold, in particular, an infinite-order homogenised 
equation is obtained by by considering a minimisation problem over a restricted set of trial fields: 


K K 


mm 


f |^^ie-^+'fc(^’^)V^(curlw-Ji)v'(curlw-Ji)| =: Ei^(w,Ji,e), 

j=0 1=0 ^ 

where is a tensor of order j + l + 2 given by 

fcOV) = (curl (curl +.^('+1))^ 

Formally considering the case K = oo leads to the infinite-order homogenised equation 

OO _ OO _ 

curl{fe(°’°)curlw(°°)} + ^e”curl{fe'‘V"curlw('^)} = ^ e^curljfe'^V^Ji}, 


( 6 . 20 ) 


n=l 


n=0 


where 




-5 E 


. -I- ... 

j+li2...ijiiij+2...ij+l+2 ' ' ' ij+li2...ijij+l+2ij+2...ij+l+iil 


j+l=n, 

iJeNufo} 


}■ 


( 6 . 21 ) 


( 6 . 22 ) 


(6.23) 


Applying the symmetrisation procedure of Section [5l it is checked that the equations (j6.17p and (j6.23p 
coincide. 

By analogy with a discussion in Section 16.31 we derive two versions of the infinite-order constitutive 
law between the electric field and electric displacement, using the Maxwell equations in combination with 
the infinite-order homogenised equation (j6.17p and by the analysis of electric energy: 


where 


j=0 n=0 j+l=n 

iJGNUfO} 


. ■= . 


(6.24) 
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Appendix A: Poincare-type inequality 

In the proofs of Theorems 12.1113.11 we use the following statement. 

Lemma 6.1. For all v € L^(T)^ such that curlv G [L2(T)]^div V G L^(T), and (v)t = 0, the inequality 

||v||[^2(t)] 3 < |T|(^||curlv||f^2(T)]3 + lldivvll^a^^^) (6.25) 

holds. 

Proof. By virtue of v G [L^(T)]" \ (v)t = 0, we write the Fourier series 

v(x) = ^ Ck G C^, k G Z^. 

kGZ3\{0} 

Noting that for k G Z^ one has curl(cke'*^'^) = i(k x Ck)e'*^'^, div(cke'*^'^) = i(k • Ck)e'*^'^, it follows that 

I|v||[l2(t)]3 = |T| < |T| 

kGZ3\{0} kGZ3\{0} 

= |T| Y1 |k-CkP) = |T|(^||curlv||f^2(T)]3 + lldivvllla^^^), 

kGZ3\{0} 

where we used the Parseval identity and the equality |kp|ckp = |k x Ckp + |k • Ckp. □ 
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